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ABSTRACT 

Aims. Intranight variation (or microvariation) is a common phenomenon of radio-loud BL Lac objects. However, it is not clear whether 
the recently found radio-quiet BL Lac objects have the same properties. The occurrence rate of intranight variation is helpful in dis¬ 
tinguishing the mechanism of the continuum of radio-quiet BL Lac objects. 

Methods. We conducted a photometric monitoring of 8 radio-quiet BL Lac objects by the Xinglong 2.16m and Lijiang 2.4m tele¬ 
scopes. The differential light curves are calculated between each target and two comparison stars. To quantify the variation, the 
significance of variation is examined by a scaled F-test. 

Results. No significant variation is found in the 11 sessions of light curves of 8 radio-quiet BL Lac objects (one galactic source is 
excluded). The lack of microvariation in radio-quiet BL Lac objects is consistent with the detection rate of microvariation in normal 
radio-quiet AGNs, but much lower than for radio-loud AGNs. This result indicates that the continua of the radio-quiet BL Lac objects 
are not dominated by jets that will induce frequent microvariations. 

Key words. Galaxies: active - Radiation mechanisms: general - BL Lacertae objects: general 


1. Introduction 

Active galactic nuclei (AGNs) are characterized by their broad 
band continua, strong emission lines, and fast variability. 
However, a handful of abnormal AGNs have recently been dis¬ 
covered in the SDSS data, i.e., weak line quasars and radio-quiet 
BL Lac objects (Diamond-Stanic et al. 2009; Plotkin et al. 2010). 
The UV/optical emission lines are absent in their UV/optical 
spectroscopies, though the shape and luminosity of their con¬ 
tinua are comparable to the normal AGNs. The fraction of these 
special AGNs is small (~1/1000) in the SDSS DR 7 sample; 
however, it could be an important stage in the evolving sequence 
of AGNs (Hryniewicz et al. 2010; Liu & Zhang 2011). In the 
early stage of an active cycle of AGNs, the radiative feedback 
can expel the gas from broad line regions and result in weak or 
even in the disappearance of broad emission lines (Liu & Zhang 
2011). Other models, such as of a cold accretion disk, an ex¬ 
tremely high accretion rate, shielding gas, or abnormal BL Lac 
objects, have also been proposed (Shemmer et al. 2010; Plotkin 
et al. 2010; Laor & Davis 2011; Wu et al. 2012). The origin of 
such weak-line AGNs is still not clear. Both the thermal (accre¬ 
tion disk) or non-thermal (jet) component may explain the weak 
line feature. 

The optical continua of classical BL Lac objects are dom¬ 
inated by the synchrotron emission from the relativistic jets, 
therefore high polarization and fast variability are important 
characteristics of classical BL Lac objects. Heidt & Nilsson 
(2011) found the polarization degrees of the radio-quiet BL Lac 
candidates are low. Plotkin et al. (2010) investigated the long¬ 
term variability of radio-quiet BL Lac objects using the data 
from SDSS stripe 82 and find that the variation amplitude of 
radio-quiet BL Lac objects is smaller than that of radio-loud BL 


Lac objects. However, owing to the small size of their sample, 
this is not very conclusive. 

The variation in short time scale (intranight) is another char¬ 
acteristic of classical BL Lac objects. During a very short pe¬ 
riod, e.g. several hours, the flux of classical BL Lac object can 
change by several tenths of a magnitude (Wagner & Witzel 1995; 
Heidt & Wagner 1996; Bai et al. 1998; de Diego et al. 1998). 
However, it is still unclear whether the intranight variation is 
frequent in the radio-quiet BL Lac objets, which will be helpful 
in distinguishing the origin of their continua. Gopal-Krishna et 
al. (2013, hereafter GJC2013) and Chand et al. (2014, hereafter 
CKG2014) claim that they have detected a considerable frac¬ 
tion of intranight variation in a sample of weak-line AGNs (duty 
cycle~5%), and this fraction can be higher if the signal-to-noise 
ratio of the light curve is further increased. However, as we dis¬ 
cuss in this paper, some galactic sources may contaminate their 
sample. 

In this paper, we report the result of our monitoring cam¬ 
paign of radio-quiet BL Lac objects. Observations and data re¬ 
duction are described in Section 2, and then the significance of 
the intranight variation is shown in Section 3. Section 4 discusses 
the implication of our results and presents the conclusions. 

2. Observations and data reduction 

The amplitude of intranight variation is normally several tenths 
of a magnitude, so the desired error of our observed mag¬ 
nitude is <0.05 mag with the exposure time not longer than 
10 min. The corresponding magnitude threshold is R ~ 18.5 
for the 2 m class telescopes we used. Seven radio-quiet BL 
Lac objects were selected from Plotkin et al. (2010), and 
SDSS J094533.99-1-100950.1 was selected from Hryniewicz et 
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Table 1. Information on observations. 


Object (SDSS) 

RA 

DEC 

R 

Redshift 

Date 

Telescope" 

Filter 

Duration (h) 

N" 

1081250.80-1-522530.8 

123.212 

52.425 

17.85 

1.152 

2011.2.12 

L 

R 

4.4 

27 






2012.1.27 

L 

R 

5.7 

32 

J085025.60-H342750.9 

132.607 

34.464 

18.51 

1.389 

2011.2.13 

L 

R 

3.3 

20 

J090107.64-H384658.8 

135.282 

38.783 

17.87 

unknown 

2012.1.29 

L 

R 

6.9 

32 

1094533.99+100950.1 

146.392 

10.164 

17.45 

1.662 

2011.2.10 

L 

V 

6.1 

25 






2012.1.28 

L 

R 

6.0 

34 

1125219.48+264053.9 

193.081 

26.682 

17.51 

1.292 

2011.4.23 

X 

R 

6.3 

38 






2013.4.14 

X 

R 

8.4 

47 

1132809.59+545452.8 

202.040 

54.915 

17.59 

2.096 

2011.4.24 

X 

R 

4.6 

48 

1134601.29+585820.1 

206.505 

58.972 

17.46 

1.667 

2011.4.25 

X 

R 

4.1 

18 

1142943.64+385932.2 

217.432 

38.992 

17.26 

0.930 

2013.4.13 

X 

R 

6.2 

33 


L-Lijiang 2.4m, X-Xinglong 2.16m 
® Number of exposures 


al. (2010). The above selection criteria are similar to those in 
GJC2013 and CKG2014. Actually, five sources in our sample 
are shared with GJC2013 and CKG2014 who based their se¬ 
lection primarily on classification by Plotkin et al. (2010) as a 
‘high-confidence BL Lac candidate’. We additionally included 
in our sample some low-confidence BL Lac candidates. The 
sources are classified as low-confidence only because the contin¬ 
uum near the emission line is hard to define, and the equivalent 
width of emission lines will be larger or smaller than 5 A de¬ 
pending on the continuum assumptions, which is mainly due to 
the noisy spectra around some emission lines. We therefore think 
there should be no systematic difference between high- and low- 
confidence sources and will investigate the variation property of 
sub samples in future works. 

The observations were carried out by BFOSC (BAG Faint 
Object Spectrograph and Camera) on the Xinglong (China) 2.16 
m telescope and YFOSC (Yunnan Faint Object Spectrograph 
and Camera ) on the Lijiang (China) 2.4 m telescope. All ob¬ 
servations were performed in Johnson R band, except for SDSS 
J094533.99-I-100950.1 in Johnson V band. The exposure time 
was 300 s or 600 s depending on the weather conditions. The 
detailed information about the sample and observations is shown 
in Table 1. In total, there are 11 sessions of light curves of these 
eight sources. 

The photometric data were reduced with the standard rou¬ 
tines in the Image Reduction and Analysis Facility (IRAF) soft¬ 
ware. The bias frames were extracted from no fewer than ten 
frames, and the flat frames did not have fewer than five frames 
in one band for one night of observation. The dark of the CCD is 
negligible (compared with the readout noise and the flat fluctua¬ 
tion) and therefore not considered. The flat frames for the same 
band were combined by average, and then the normalized flat 
frame was generated; the normalized bias frame was generated 
by median combination. Then the source frames were corrected 
by the normalized bias frame and flat frame. 

With the corrected source images, we used the package 
APPHOT to perform aperture photometry. The values of en¬ 
closed, moffat, and direct for the comparison stars and the target 
source were used to estimate the mean full width at half maxi¬ 
mum (FWHM). The apertures of the photometry for individual 
frame were carried with 2.5 ~ 3 times of FWHM. If the value of 
FWHM significantly changed during one night, we took differ¬ 
ent values of FWHM even for the same source. 


3. Results 

To detect the underlying variation of the target, we first cal¬ 
culated the differential light curves (DLCs) between the target 
and comparison stars. Two nearby comparison stars (noted as 
Star 1 and Star 2 hereafter) with similar magnitudes to the target 
were selected and the DLCs of AGN-Star 1, AGN-Star 2, and 
Star 1 -Star 2 are shown in Figure 1. The position and g-r color 
of targets and comparison stars are shown in Table 2. Since the 
color difference between target and star pairs is smaller than 1.5, 
the variation in air mass during the observation has little effect on 
DLCs (Carini et al. 1992; Stalin et al. 2004). We also tried dif¬ 
ferent companion stars, and the final significance of intranight 
variation is quite robust. Some exposures with bad weather were 
excluded from the DLCs, which led to some gaps in the DLCs. 

To quantify the significance of the variation of light curves, 
we performed a scaled F-test, which is more powerful and re¬ 
liable than the traditional C-test (de Diego 2010). The scaled F 
value (Howell et al. 1988) is defined as 

77 ■'AGN-Star 1 /,n 

^ = 

^ ‘^Star 1-Star 2 
N _ 2 

where ^ ^ can stand for AGN-Star 1 or 

i-l 

Star 1-Star 2. 

The definition of F^ is 


/ A^Star 2 

■^Star l^-^AGN + ^) + A^AGN^-^Star 1 

\ A^agn / 

2(^S.ar 1 + P) + l(^Star 2+P)\ 


which is the scaled factor to account for the different accura¬ 
cies between the photometries of the target and comparison stars 
(Howell et al. 1988). The variables Nagn, A^star i, and Astar 2 are 
the total counts (sky-subtracted) of target, Star 1, and Star 2, re¬ 
spectively. The variable P is defined as P = np(Ns + N}), where 
Hp is the number of pixels in the applied measuring aperture, the 
variable As is the sky photons per pixel, and Ny is the readout 
noise (e“/pixel). The value of F^ can be calculated frame-by- 
frame. However, the variation in F^ of our observations during 
one night is no more than 10% owing to the small variation of 
our targets. Therefore, we have taken the median value of F^ for 
the exposures in one night. Our final result is not sensitive to this 
choice. 

The significance of the variation is determined by the F dis¬ 
tribution with AAGN-Star 1 - 1 and Astar I-Star 2 “ 1 degrees of 
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SDSS J081250.80+522530.8 (2011.2.12) SDSS J081250.80+522530.8 (2012.1.27) 




SDSS J085025.60+342750.9 (2011.2.13) 



SDSS J090107.64+384658.8 (2012.1.29) 



14 15 16 17 18 19 20 21 


UT (h) 


SDSS J094533.99+100950.1 (2011.2.10) SDSS J094533.99+100950.1 (2012.1.28) 




Fig. 1. Differential light curves of AGN-Star 1, AGN-Star 2, and Star 1-Star 2. 
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SDSS J125219.48+264053.9 (2011.4.23) 



14 15 16 17 18 19 20 

UT (h) 


SDSS J125219.48+264053.9 (2013.4.14) 



SDSS J132809.59+545452.8 (2011.4.24) SDSS J134601.29+585820.1 (2011.4.25) 




SDSS J142943.64+385932.2 (2013.4.13) 
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Table 2. Information on targets and their companion stars. 


Object (SDSS) 

Date 

RA (J2000) 

DEC (J2000) 

r (SDSS) 

g - r (SDSS) 

J081250.80+522530.8 

2011.2.12 

08 12 50.80 

+52 25 30.8 

18.05 

0.3 

Star 1 


08 12 51.29 

+52 26 46.4 

17.28 

1.4 

Star 2 


08 12 49.52 

+52 26 26.2 

17.89 

1.3 

1081250.80+522530.8 

2012.1.27 

08 12 50.80 

+52 25 30.8 

18.05 

0.3 

Star 1 


08 12 51.29 

+52 26 46.4 

17.28 

1.4 

Star 2 


08 12 49.52 

+52 26 26.2 

17.89 

1.3 

J085025.60+342750.9 

2011.2.13 

08 50 25.60 

+34 27 50.9 

18.66 

0.4 

Star 1 


08 50 26.96 

+34 26 35.9 

19.22 

1.0 

Star 2 


08 50 17.77 

+34 26 50.5 

18.53 

0.7 

J090107.64+384658.8 

2012.1.29 

09 01 07.64 

+38 46 58.8 

18.12 

0.1 

Starl 


09 01 06.48 

+38 47 08.7 

18.14 

1.4 

Star 2 


09 01 05.15 

+38 48 24.5 

17.36 

1.3 

1094533.99+100950.1 

2011.2.10 

09 45 33.99 

+ 10 09 50.1 

17.66 

0.4 

Star 1 


09 45 27.96 

+ 10 08 47.8 

16.89 

0.4 

Star 2 


09 45 37.93 

+ 10 08 08.9 

18.01 

0.7 

1094533.99+100950.1 

2012.1.28 

09 45 33.99 

+ 10 09 50.1 

17.66 

0.4 

Star 1 


09 45 27.96 

+ 10 08 47.8 

16.89 

0.4 

Star 2 


09 45 37.93 

+ 10 08 08.9 

18.01 

0.7 

J125219.48+264053.9 

2011.4.23 

12 52 19.48 

+26 40 53.9 

17.70 

0.2 

Starl 


12 52 27.12 

+26 38 49.7 

17.44 

1.1 

Star 2 


12 52 14.26 

+26 39 11.5 

17.15 

1.3 

J125219.48+264053.9 

2013.4.14 

12 52 19.48 

+26 40 53.9 

17.70 

0.2 

Starl 


12 52 23.02 

+26 38 42.9 

15.82 

0.6 

Star 2 


12 52 23.82 

+26 41 42.6 

16.43 

0.3 

J132809.59+545452.8 

2011.4.24 

13 28 09.59 

+54 54 52.8 

17.84 

0.1 

Star 1 


13 27 58.21 

+54 54 00.2 

17.54 

0.8 

Star 2 


13 28 22.83 

+54 55 54.7 

18.14 

0.5 

1134601.29+585820.1 

2011.4.25 

13 46 01.29 

+58 58 20.1 

17.74 

0.2 

Starl 


13 46 06.60 

+58 58 08.2 

18.01 

1.5 

Star 2 


13 45 55.76 

+58 57 34.8 

18.46 

1.3 

J142943.64+385932.2 

2013.4.13 

14 29 43.64 

+38 59 32.2 

17.55 

0.0 

Starl 


14 29 39.99 

+39 02 19.6 

17.21 

0.9 

Star 2 


14 29 30.47 

+39 00 08.7 

16.33 

0.6 


freedom, where A^AGN-star i and A^star i-star 2 are the number of 
observations in the AGN-Star 1 and Star 1-Star 2 DLCs, respec¬ 
tively. 

The results of the significance are listed in Table 

3. Since we exchanged the position of Star 1 and 
Star 2 in equations (1) and (2), there are two values 
of significance for (AGN - Star 1)/(Star 1 - Star 2) and 
(AGN - Star 2)/(Star 2 - Star 1). 

As indicated by the results of T-test, we only detect a sig¬ 
nificant variation (~3cr level) in SDSS J090107.64H-384658.8. 
However, due to the large proper motion of this source (62+11 
mas/yr from Monet et al. 2003), its extragalactic nature is 
doubtable. Therefore, we would like to exclude it from the fi¬ 
nal sample of radio-quiet BL Lac objects. As a result, there is no 
significant variation detected in our observations of radio-quiet 
BL Lac objects. 

4. Discussions and conciusions 

Radio-loud AGNs and blazars can exhibit microvariation with 
a large amplitude up to ~100%. However, some microvariation 
events are also observed in radio-quiet AGNs with high signifi¬ 
cance (Stalin et al. 2004; Gupta & Joshi 2005). The mechanism 
of the microvariation in radio-loud AGNs is believed to be the 
fluctuation caused by the shocks in jets. However, the instabil¬ 
ity or flares in the accretion disk can also induce microvariation 
even for the radio-quiet AGNs (Mangalam & Wiita 1993). A 


weak blazar component in radio-quiet AGNs is an alternative to 
microvariation (Czerny et al. 2008). Though the occurrence of 
microvariations is not a smoking gun of jets, the fraction and 
amplitude of the microvariations in radio-quiet and radio-loud 
ones are quite different. 

Gupta & Joshi (2005) compiled the microvariations of dif¬ 
ferent classes of AGNs and found the detection fractions of mi¬ 
crovariation in radio-quiet and radio-loud (non-blazars) AGNs 
are ~10% and ~35-40%, respectively. For blazars, the fractions 
are ~60-65% and ~80-85% for the observations that are less than 
and more than 6 h, respectively. In addition, they also claim that 
the amplitude of the microvariation of radio-loud ones is larger 
than that of radio-quiet ones. 

Carini et al. (2007) established a sample of 117 radio-quiet 
AGNs that have been investigated for microvariations and found 
a detection rate of microvariations for the entire sample of 
21.4%. If the criteria for ‘radio-quiet’ are strengthened to R < 1 
(R is the ratio of the radio [5 GHz] flux to optical [4400 A] flux), 
the detection rate of microvariations is only 15.9%. 

Goyal et al. (2013) analyzed 262 sessions of light curves 
of 77 AGNs from their uniform AGN monitoring data and 
found the duty cycles of intranight variation of radio-quiet 
quasars, radio-intermediate quasars, lobe-dominated quasars, 
low optical-polarization core-dominated quasars, high optical- 
polarization core-dominated quasars, and TeV blazars are 10%, 
18%, 5%, 17%, 43%, and 45%, respectively. 
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Table 3. Results of the significance of variations. 


Object (SDSS) 

Date 

Ti 

Tj 

Li 

F2 

Significance! 

Significance 2 

1081250.80+522530.8 

2011.2.12 

1.16 

1.66 

0.70 

1.34 

18.4% 

76.9% 

1081250.80+522530.8 

2012.1.27 

1.50 

1.85 

0.51 

0.72 

3.44% 

18.7% 

1085025.60+342750.9 

2011.2.13 

1.09 

0.59 

1.20 

1.41 

64.9% 

77.0% 

1090107.64+384658.8 

2012.1.29 

1.84 

1.34 

1.75 

2.68 

93.8% 

99.6% 

1094533.99+100950.1 

2011.2.10 

0.45 

1.25 

1.15 

0.93 

63.5% 

43.0% 

1094533.99+100950.1 

2012.1.28 

0.51 

1.26 

0.74 

1.16 

20.1% 

66.2% 

1125219.48+264053.9 

2011.4.23 

1.80 

1.58 

0.39 

0.36 

0.25% 

0.13% 

1125219.48+264053.9 

2013.4.14 

6.76 

7.15 

0.49 

0.46 

0.91% 

0.50% 

1132809.59+545452.8 

2011.4.24 

0.79 

1.29 

1.05 

0.91 

56.6% 

37.3% 

1134601.29+585820.1 

2011.4.25 

0.60 

0.96 

0.67 

1.18 

20.9% 

63.1% 

1142943.64+385932.2 

2013.4.13 

2.05 

1.40 

1.51 

2.23 

87.4% 

98.7% 


Notes. The subscripts ‘L and ‘2’ of variables stand for the results of (AGN - Star l)/(Star 1 - Star 2) 
and (AGN - Star 2)/(Star 2 - Star 1), respectively. 


No significant microvariation is detected in our final sample, 
in ten sessions of light curves. The Icr upper limit of the fraction 
of microvariation is 15% using the method of Cameron (201 ifl 
In deriving this upper limit, we treated the sources equally. The 
weights of sources should not be the same owing to different 
exposure times, signal-to-noise ratios, and observation numbers; 
however, this potential minor correction will not change our fi¬ 
nal conclusion. This low fraction in our sample of radio-quiet 
BL Lac objects is consistent with that of the radio-quiet AGNs 
but much lower than for the radio-loud ones and blazars. This 
indicates that the continuum of radio-quiet BL Lac objects is not 
dominated by the jet component. Actually, the SED of radio¬ 
quiet BL Lac objects is similar to the normal radio-quiet AGNs 
(Lane et al. 2011), which further supports the accretion disk ori¬ 
gin of the continuum. Accurate black hole mass measurements 
can determine the accretion state of radio-quiet BL Lac objects 
and further distinguish the different models related to the accre¬ 
tion disk origin of their continua, which will be explored in our 
future works. 

GJC2013 and CKG2014 detected signihcant variations 
(confidence level >99% for two comparison stars) of SDSS 
J090843.25-H285229.8 and SDSS J121929.45-I-471522.8 in their 
29 light curves. Based on these two events, they derived a 
duty cycle ~5% of intranight variation from their sample on 
weak-line AGNs. However, two sources in their sample (SDSS 
J090107.64-H384658.8 and SDSS J121929.45 h- 471522.8) are 
likely to be galactic sources owing to the large proper mo¬ 
tion. The proper motions of SDSS J090107.64-1-384658.8 and 
SDSS J121929.45-I-471522.8 from USNO-B are 62+11 mas/yr 
and 112+4 mas/yr, respectively (Monet et al. 2003). These two 
bright sources (y~18) are well above the completeness limit 
y=21 of the USNO-B catalog (Monet et al. 2003). The posi¬ 
tional error of one epoch is ~200 mas. Thus, the proper motion 
of ~100 mas/yr can be accurately detected with the epoch dif¬ 
ference of ~40 years. The systematic error of the proper mo¬ 
tion is comparable to the statistical error. Munn et al. (2004) and 
Roeser et al. (2010) have further calibrated the USNO-B cata¬ 
log with SDSS and 2MASS astrometry, and the resulting proper 
motions of these two sources are consistent with those from the 
USNO-B catalog. No radio or X-ray counterpart is found near 
their positions. Actually, Rebassa-Mansergas et al. (2010) lists 

' Given the sample size n and observed success counts k, the upper 
limit pu is defined by (a- 1 )*(/,‘- = (1 “ 0/2, where a = 
k+\,b = n- k+ \,q = 1-p, and c is the confidence level. 


SDSS J121929.45+471522.8 as a DC white dwarf in their cat¬ 
alog. SDSS J090107.64+384658.8 is classified as an uncertain 
DC white dwarf by Eisenstein et al. (2006) and as an uncertain 
DC+M binary system by Kleinman et al. (2013). The variation 
observed is likely due to the oscillation or accretion of the white 
dwarfs (Winget & Kepler 2008; Eontaine & Brassard 2008). 

The DC white dwarfs are generally difficult to positively 
identify due to their featureless spectra and the discrepancy in¬ 
deed exists in different catalogs. However, we think these pos¬ 
sible galactic sources should be excluded from the sample of 
AGNs to be safe. We detected the intranight variation in SDSS 
J090107.64+384658.8 but excluded it from the sample. After 
the possible white dwarfs are removed from the sample of 
GJC2013 and CKG2014, only one significant variation (SDSS 
J090843.25+285229.8) is detected in the remaining 22 light 
curves, which is also the only significant event in all 32 light 
curves (including 10 additional ones of this paper) of weak-line 
AGNs up to now. This low occurrence rate is consistent with the 
rate for the normal radio-quiet AGNs. 

Our present sample of radio-quiet BL Lac objects is still 
too small to constrain the duty cycle of the microvariation well. 
We will enlarge our sample, especially for the monitoring time 
longer than six hours, and improve the accuracy of photometry 
to ~ 0.01 mag to detect smaller variations. 

Observations in more bands will help to investigate the color- 
behavior and further constrain the mechanism of the continuum 
of the radio-quiet BL Lac objects. 
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